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OPTICAL LINE DIAGNOSTICS OF Z «2 OPTICALLY FAINT ULTRA-LUMINOUS INFRARED GALAXIES IN THE 
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Abstract 

We present near-infrared spectroscopic observations for a sample of ten optically faint luminous infrared 
galaxies (R- [24] > 14) using Keck NIRSPEC and Gemini NIRI. The sample is selected from a 24 /im Spitzer 
MIPS imaging survey of the NDWFS Bootes field. We measure accurate redshifts in the range 1.3 <z<3.4. 
Based on either emission line widths or line diagnostics, we find that all ten galaxies harbor luminous AGN. 
Seven sources are type I AGN, exhibiting broad (> 1900 km s" 1 ) Ha or H/3 emission lines; the remaining three 
are type II AGN. Given their large mid-IR luminosities and faint optical magnitudes, we might expect these 
sources to be heavily extincted quasars, and therefore only visible as type II AGN. The visibility of broad lines 
in 70% of the sources suggests that it is unlikely that these AGN are being viewed through the mid-plane of a 
dusty torus. For four of the sources we constrain the Ha/H/3 Balmer decrement and estimate the extinction to 
the emission line region to be large for both type I and type II AGN, with A^a ^2.4-5 mag. Since the narrow- 
line region is also extincted and the UV continuum emission from the host galaxies is extremely faint, this 
suggests that much of the obscuration is contributed by dust on large (~kpc) scales within the host galaxies. 
These sources may be examples of "host-obscured" AGN which could have space densities comparable or 
greater to that of optically luminous type I AGN with similar bolometric luminosities. 

Subject headings: galaxies: active — galaxies: starburst — infrared: galaxies — quasars: general — quasars: 
emission lines 



1. INTRODUCTION 

Recent results from the Spitzer Space Telescope (Werner 
et al. 2004) have shown that luminous infrared galaxies 
(LIRGs) become increasingly important at high redshift (Le 
Floc'h et al. 2004). They dominate the energy density of the 
Universe at z > 1 and their extrapolated infrared luminosi- 
ties suggest that they are galaxies undergoing a highly active 
phase in the growth of their bulges and / or accompanying 
super-massive black holes (SMBHs). Although their power 
source appears to be primarily due to starbursts at lower lu- 
minosities (Li R (8- 1000 pm) < 10 12 L Q ), AGN activity be- 
comes increasingly important at higher luminosities (Veilleux 
et al. 1995; Dudley 1999; Brand et al. 2006b). 

We have uncovered a population of sources which are bright 
at 24 pm but optically very faint (Dey et al. in prep.). Al- 
though only a relatively small fraction of the total 24 pm pop- 
ulation, the fraction increases from 5% at 24 /im flux densi- 
ties (f24) greater than 1 mJy to 15% at f^n > 0.3 mJy (Dey 
et al. in preparation), suggesting that they become increas- 
ingly important either at high redshifts or low luminosities. 
Hereafter, we will refer to sources with R— [24] > 14 as op- 
tically faint luminous infrared galaxies. Optically faint lu- 
minous infrared galaxies have optical to mid-infrared colors 
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which are more than a magnitude redder than Arp 220 (an 
extreme dusty source in the local Universe). Their huge bolo- 
metric luminosities imply that they must be undergoing rapid 
growth via dust-enshrouded AGN and/or starburst activity. 
However, because these sources are also optically faint, they 
are also likely very dust-obscured and will have been largely 
missed by previous optical surveys. Determining the evolu- 
tion of their luminosity function and the relative contribution 
of AGN and starburst activity to their bolometric luminosity 
is critical in understanding how the most massive galaxies and 
their SMBHs built up. 

Although optically faint luminous infrared galaxies are dif- 
ficult to study at optical wavelengths, their extreme red SEDs 
makes them good candidates for mid-infrared (mid-IR) spec- 
troscopy. Our follow-up of a sub-sample of 58 of these 
sources with Spitzer I IRS, yielded redshifts in the range 
z w 2 - 3 for 34 of the sources, based primarily on the strong 
9.7 /jm silicate absorption feature (Houck et al. 2005; Weed- 
man et al. 2006; Higdon et al. in preparation; see also Yan 
et al. 2005). An additional 9 have probable redshifts derived 
from weak silicate absorption features (Weedman et al. 2006). 
15 sources have featureless power-law IRS spectra (Weedman 
et al. 2006) from which one can obtain little further informa- 
tion. It is not clear whether this is due to the sources lying 
at redshifts for which the silicate absorption feature falls out 
of the observable band (z > 2.5), or because they have only 
weak silicate absorption features (Weedman et al. 2006). In 
these extreme sources, the IRS spectra often have low signal- 
to-noise ratios and can be difficult to interpret. Usually, the 
silicate absorption feature is the only strong observable fea- 
ture. The strength of this feature is determined primarily by 
the geometry (Levenson et al. 2006) or dumpiness (Spoon 
et al. 2006) of the dust region and can be strong in both AGN- 
and starburst-dominated sources (although Spoon et al. 2006 
show that starburst-dominated sources tend to have stronger 
PAH features and only shallow absorption features). The sili- 
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cate strength may simply provide us with a lower limit to the 
line of sight obscuration towards the nucleus. Redshifts and 
diagnostics from other wavelengths will help us understand 
these sources better and may be able to identify starburst sig- 
natures that may be overwhelmed by the hot dust continuum 
in the mid-IR. 

Rest-frame optical spectra can potentially help in under- 
standing the nature of these sources. One can determine red- 
shifts for the sources with featureless power-law spectra in 
the mid-IR, and confirm and refine the redshift measurements 
for those sources with broad silicate absorption features. By 
determining the widths and relative strengths of the optical 
lines, one can also help constrain the properties of the energy 
source. Ha is a particularly useful line in this regard because 
it is generally luminous and it can be significantly broadened 
by the presence of an AGN. In the cases in which the broad- 
line region is obscured from view, one can use the line ratios 
of [O III] A5007 / H/3 and [Nil] A6583 / Ha to distinguish 
between narrow-line AGN and H II region-like galaxies (e.g., 
Baldwin et al. 1981; Veilleux & Osterbrock 1987; Veilleux 
et al. 1995; Kauffmann et al. 2003; Kewley et al. 2006). The 
/?-[24]>14 selection criteria means our sample should be par- 
ticularly affected by dust extinction. Ha should be less at- 
tenuated than emission lines at shorter wavelengths and the 
Balmer decrement (i.e., the Ha to H/3 ratio) provides a useful 
measure of the optical extinction. 

For the redshifts of the sources in our sample, the typically 
strong, diagnostically useful optical lines of Ha and H/3 are 
shifted into the near-IR. Being very red, optically faint lu- 
minous infrared galaxies can also be relatively bright in the 
near-IR (K w 19-20) compared to the optical (R > 24). Near- 
IR spectroscopy using the largest telescopes presents the best 
(and in many cases, the only) chance of determining the na- 
ture of these extreme sources. In this paper, we present and 
analyze Keck NIRSPEC and Gemini NIRI near-IR spectra of 
10 high redshift optically faint ULIRGs. We use optical emis- 
sion line diagnostics to determine their redshifts and provide 
insights into their primary power source. 

The paper is structured in the following way. The sample 
selection criteria are presented in Section 2. In Section 3, we 
discuss the observational techniques, data reduction, and cal- 
ibration. The near-IR spectra are presented in Section 4. We 
also discuss each source individually. In Section 5, we discuss 
the characteristics of the sample. Section 6 is a discussion sec- 
tion in which we investigate the broader implications of our 
results. A cosmology of Hq = 70 km s" 1 Mpc" 1 , Q M =0.3, and 
Sl\=0.7 is assumed throughout. 

2. THE SAMPLE 

Our sample of 10 optically faint ULIRGs is selected from a 
24 /im survey of the NOAO Deep Wide-Field Survey (ND- 
WFS; Jannuzi & Dey 1999) Bootes field. The NDWFS 
Bootes field comprises deep optical (/?<25.5) and near-IR 
(K < 18.6) imaging. A thorough description of the NDWFS 
survey will be provided in Jannuzi et al. and Dey et al. (both 
in preparation). The ^9.3 deg 2 region has been mapped by 
Spitzer MIPS at 24 /im (f^) and comprises «20,000 sources 
down to a 5 (J depth of «0.3 mJy (Le Floc'h et al. in prepa- 
ration). For more details regarding the creation of the 24/im 
catalog, see Brown et al. (2006) and Brand et al. (2006b). The 
^2,500 R- [24] > 14 sources will be described in more detail 
in Dey et al. (in preparation). 

Our sample is selected from the R- [24] > 14 sources using 
one of the following three criteria: 



(1) Sources with a known redshift in the range 1.25 < z < 
1.6, 2.1 < z < 2.6, or 3.2 < z < 3.7. At these red- 
shifts, the strongest expected emission lines (i.e., Ha 
/ [Nil] and H/3 / [OIII] ) have rest-frame wavelengths 
which are observable from the ground. Although this 
sample is biased to sources for which we could ob- 
tain redshifts (i.e., sources with strong UV emission 
lines observable with optical spectroscopy and/or a 
strong silicate absorption feature within the Spitzer IRS 
wave-band), it provided us with objects for which the 
near-IR spectra would yield detections or useful con- 
straints on the rest-frame optical lines. This sample 
consists of 2 sources for which we obtained optical red- 
shifts (Desai et al. in preparation) but which exhibited 
featureless power-law spectra in the mid-IR (SST24 
j 143424.4+334543 and SST24 J142644.3+333051), 1 
source with an optical redshift with no IRS observation 
(SST24 J14301 1.3+343439), and 3 sources for which 
the silicate feature in the IRS spectra provided an esti- 
mate of the redshift (SST24 J 142827. 1+354127, SST24 
J143312.7+342011,andSST24J143028.5+343221). 

(2) Sources observed by Spitzer IRS which yielded only 
featureless mid-IR spectra and no redshift measure- 
ment. Observing these sources with near-IR spec- 
troscopy may be the only way to obtain a redshift and 
determine their nature. Four such sources were tar- 
geted by our near-IR spectroscopy, but only two had 
detected emission lines and are presented here (SST24 
J142842.9+342409 and SST24 J142939. 1+353558). 

(3) Sources which are bright in the /f-band (K w 17.5- 
18.0). Six such sources were selected and tar- 
geted by our Gemini NIRI observations but only two 
sources exhibited strong emission lines and are pre- 
sented here (SST24 J143027.1+344007 and SST24 
J 142800.7+350455). 

In all cases without a bright K-bwA detection, we preferen- 
tially chose sources that had a K < 18 star within 20". This 
made target acquisition quicker and easier and provided a ref- 
erence from which to measure where spectrum of the target 
should fall (see section 3). We also chose sources with the 
brightest detections in the TT-band. Although the TT-band light 
may not be correlated with the line emission strength, detect- 
ing continuum is useful for confirming that the source is in the 
slit and for providing a relative flux calibration of the source. 
The /f-band magnitudes were obtained from either our ND- 
WFS ONIS observations, the FLAMEX survey (Elston et al. 
2006), or follow-up Keck NIRC observations (Weedman et al. 
2006). 

A table summarizing our sample and their auxiliary data is 
provided in Table 1 . 

3. OBSERVATIONS AND DATA REDUCTION 

Our sample was observed with either the Near Infrared 
Spectrograph (NIRSPEC; McLean et al. 1998) on the Keck 
II telescope or the Near Infrared Imager (NIRI; Hodapp et al. 
2003) on the Gemini North telescope. We discuss the observa- 
tions separately in the following two sections and then discuss 
the data reduction which was similar for both data sets. 

3.1. Keck NIRSPEC 

The majority (8/10) of our sample were observed with NIR- 
SPEC on the Keck II telescope on 26th - 27th May 2005 
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TABLE 1 

Observed and derived properties of the sample of 10 optically faint ULIRGs. 



Name" J142827.1+354127 J143027. 1+344007 J1433 12.7+342011 J14301 1.3+343439 b J143028.5+343221 



RA (J2000) c 


14:28:27.191 


14:30:27.203 


14:33:12.708 


14:30:11.367 


14:30:28.521 


Dec (J2000) 


35:41:27.70 


34:40:07.86 


34:20:11.20 


34:34:39.68 


34:32:21.34 


B w mag (Vega) d 


23.85 


25.63 


24.62 


22.88 


25.17 


R mag (Vega) d 


22.78 


24.76 


24.37 


22.02 


24.45 


/ mag (Vega) d 


21.64 


23.19 


23.57 


21.84 


24.16 


f24 M m (mJy) 


10.55 


1.17 


1.76 


1.12 


1.27 


R-[24] 


15.7 


15.3 


15.3 


12.5 


15.1 


1Kb z 


i n 
1.33 




L.L 




L.L 


Optical z c ^ 


i.zyz 






I. YL 


2.1 /O 


Near-IR instrument / band 




MTPT / T/TJ 
1N1K1 / J 111 








Exp (s) 


JOUU 


JJDU / Jo^fU 


JDUU 


3 OUU/ J OUU 


JOUU 




9ft05mav9fi 


?nf)fiinnnR/07 
i.yj\)\i \ \xn\joi\j i 


?006mavl 1 






Near-IR z 


1.292+0.001 / 1.383±0.001 


1.370+0.001 


2.114+0.001 


2 1 1 4+0 00 1 


2 178+0 001 


Continuum amplitude 8 


1.10+0.00 


0.59+0.01 


0.21+0.04 


0.23+0.05 


0.36+0.05 


Ha A (A) 


15045 / 15689 


15557 


20434 


20436 


20857 


Ha flux h 


19.5±5.3/115.7±12.3 


390.7+12.8 


375.1+51.8 


101.3+13.1 


104.5+18.3 


Ha <r n (A) 

L(Ha) k (10 42 erg s -1 ) 


4.8±0.9/ 86.8+5.2 


42.8+0.7 


90.8+8.3 


13.6+1.1 


16.0+1.5 


0.19+0.04/1.3+0.1 


4.4+0.1 


12.4+1.7 


3.4+0.4 


3.7+0.6 


FWHMCHa^ (km s"h 


227+42/3919+233 


1940+34 


3138+288 

Jl JO _LfcuO 


469+37 


542+49 




7 7-L2 1 /44 1 +4 7 

/ . / _1_Z.. l/'-r'T. 1 _1_*t. / 


280+10 


3 OJ _l_ I -JVJ 


142+38 


90+20 


TNTT1 A'(A^ 


1 5092 / 1 5736 


15606 


20499 


20501 


20922 


TNTT1 flux h 


< 98 8/< 64 8 


< 56.2 


5.2d=13.3 


33 2+17 6 




[1N11J a (ft.) 


— /— 






1JJIJJ 


14 4±2 2 


risjni / 
[iniij / n i v 


<^ J. 1/^*^ yJ.J 1 ) 


y*^ V. 1*+} 




n 3*1-1-0 1 s 

U.JJXU. l o 




Continuum amplitude^ 




n iQ-t-n ni 








H/3 A (A) 




1 1521 




15138 




H/3 flux h 




< 17.6 


— 


< 11.9 


— 


H/3 o- n (A) 


— 


— 


— 


— 


— 


L(H/3) k (10 42 erg s -1 ) 


— 


— 


— 


— 


— 


FWHM(H/3) (kms^ 1 ) 


- 


- 


- 


- 


- 


EW(H/3) (rest) 












[OIII] A" (A) 




11866 




15591 




[OIII] flux h 




2.8+2.8 




37.6+10.8 




[OIII] cr n (A) 




6.3+3.7 




22.5+3.5 




L([OIII]) k (10 42 ergs-') 




0.03+0.03 




1.2+0.3 




FWHM([OIII]) (kms" 1 ) 




375+220 




(1020+160) 




[OIII] / H/3 




(> 0.6) 




> 3.1 




Ha /H/3 




>22.5 




> 8.5 




E(B-V) 




>1.93 




> 1.00 




A(Ha) 




>4.64 




> 2.40 




class 1 


SB/AGN1 


AGN1 


AGN1 


AGN2 


AGN2 


L IR (lO 13 L ) m 


2.8 


0.3 


1.3 


0.8 


1.0 


IRS spectrum 


Desai et al. 2006 




Higdon et al. in prep. 




Higdon et al. in prep. 


IRS features 


Si Abs 




Si Abs 




Deep Si Abs 



UT and 10th - 11th May 2006 UT in photometric condi- 
tions. NIRSPEC is designed to operate over the wavelength 
region 0.94 to 5.4 /tm. We observed with both the TT-band 
(NIRSPEC 7) and H-band (NIRSPEC 5) filters. The observ- 
ing band and date for each individual object is given in Ta- 
ble 1. Observations were obtained in low-resolution mode 
(R«1500) with a 0"76 wide and 42" long slit. Typical seeing 
was 0"7 on the 26th - 27th May 2005 UT and 0"5 on the 10th 
-11th May 2006 UT. 

Our sources are typically too faint to acquire with short in- 
tegrations on the acquisition camera, SCAM. Instead, we took 
2 minute exposures with the K' filter on SCAM, aligned the 
slit on a nearby (<20"away) brighter source and performed a 
blind offset to the target. In a few cases in which the optical 
images showed the target to be elongated, we used a specific 
position angle to obtain the maximum amount of light down 
the slit. In all other cases, we set the position angle such that 
the offset star was also in the slit. In cases in which the target 



was very faint, this helped to determine exactly where the tar- 
get should fall on the array. For sky subtraction, the observa- 
tions were split into 10 minute exposures, and the object was 
nodded along the slit by 10"-20"(taking care that the position 
of the brighter object used for acquisition did not coincide 
with any of the dither positions of the target). 

Four sources which were known to have 2.1 < z < 2.4 from 
either the optical or IRS spectra were observed in the TT-band 
to obtain measurements of the Ha and [Nil] emission lines. 
Four sources with no redshift information were also observed 
in the /T-band (we only present the results of the 2 sources 
which were detected; the remaining 2 sources are discussed 
in Section 5.4). For 2 sources exhibiting strong Ha emis- 
sion, we also observed the sources in the //-band to deter- 
mine line strengths and widths of the [OIII] AA5007, 4959 
doublet and H/3 lines. SST24 J142827. 1+354127 is known to 
he at z =1.292 and was observed in the //-band to target Ha 
and [Nil]. SST24 J142644.3+333051, known to havez=3.355 
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TABLE 1 
Continued. 



Name a 


J142842.9+342409 


J142800.7+350455 


J143424.4+334543 


J142939. 1+353558 


J142644.3+333051 


RA (J2000) c 


14:28:42.913 


14:28:00.724 


14:34:24.500 


14:29:39.159 


14:26:44.330 


Dec (JzUUU) 


i /i ■ 9 /i ■ no /i/i 
J4.z4.Uy. 44 


9J.U4.JJ. 31 


jj.4j.4j.4U 


Q^-Q^.^Q 9/1 
Jj.Jj.j5.z4 


J J. JU. JZ.UO 


By/ mag (Vega)^ 


93 ^7 
ZJ.J / 


9s nn 


(Z4.J/J 


/97 Q7~i 

(zi.il) 




mag (Vega)^ 


99 

ZZ. J4 


94 98 
Z4.ZO 


9S £9 
ZJ.OZ 




94 T.A 
Z4. J4 


/ mag (Vega)* 1 


21 47 


91 SQ 
z J.ov 


^Z+. 3\J) 


/9S 7CA 
^ZJ. la) 






3.11 


0.92 


0.86 


1.16 


1.14 


K-[24] 


14.1 


14.5 


15.6 


17.4 


14.8 


IRSz 6 


? 




? 


? 


? 


Optical z e - f 


9 


- 


2.26 


? 


3.355 


Near-IR instrument / band 


NIRSPEC / KIH 


NIRI / HIK 


NIRSPEC / K 


NIRSPEC / K 


NIRSPEC / K 


Exp (s) 


3600 / 7200 


3840 / 960 


3600 


3600 


3600 


Date (UT) 


2006may 10/11 


2006jun08/jul29 


2006mayll 


2006mayll 


2006mayl0 


Near-IR z 


2.180±0.0()1 


2.223±0.001 


2.263±0.001 


2.498±0.001 


3.312±0.001 


Continuum amplitude 6 


0.32±0.07 


0.00±0.01 


0.30±0.03 


0.10±0.03 


- 


Ha A (A) 


20870 


21152 


21411 


22957 


- 


Ha flux h 


673.2±41.9 


230.3±27.8 


27.0±11.8 


200.0±31.1 


- 


Ho cr n (A) 


93.6±3.2 


114.9±6.9 


21.0±4.6 


67.1±5.7 


- 


L(Ha) k (10 42 ergs"') 


24.0±1.5 


8.6±1.0 


1.1 ±0.5 


9.9±1.5 


- 


FWHM(Ha)(kms _1 ) 


3167±110 


3833±229 


692±153 


2063±176 


- 


EW(Ho) (rest) (A) 


672±154 


- 


27.7±12.5 


564±189 


- 


[NII]A'(A) 


20936 


21219 


21479 


23029 


- 


[Nil] flux h 


< 103.7 


< 62.1 


124.9±16.1 


<158.3 


- 


[Nil] cr"(A) 


- 


- 


31.9±2.6 


- 


- 


[Nil] / Ha 


(< 0.15) 


« 0.3) 


4.6±2.1 


« 0.79) 


- 


Continuum amplitude 6 


0.01±0.01 


O.lliO.01 


- 


- 


1.07±0.02 


H/3 A (A) 


15459 


15669 


- 


— 


20964 


H/3 flux h 


< 42.9 


< 59.8 


_ 


_ 


57.9±17.7 


H/3 <r n (A) 

L(H/3) k (l() 42 erg s~') 










63.4±11.8 










5.7±1.7 


r vv nivi^n./j y ^Jvin » } 






















12.5±3.8 


[OIII] AJ(A) 
[OIII] flux h 


15922 


16138 






21591 


104.8±3.6 


88.2±6.() 






108.2±10.5 


[OIII] o-"(A) 


36.5±0.7 


13.4±0.4 






16.3±0.8 


L([OIII]) k (10 42 ergs-') 


3.7±0.1 


3.3±0.2 






10.6±1.0 


FWHM([OIII]) (kms- 1 ) 


1617±31 


585±19 






533±27 


[OIII] / H/3 


(> 2.44) 


(>1.5) 






(1.87±0.6) 


Ha /H/3 


>15.7 


>3.8 








E(B-V) 


>1.59 


>0.22 








A(Ha) 


> 3.81 


> 0.54 








class 1 


AGNl 


AGNl 


AGN2 


AGNl 


AGNl 


L IR (lO 13 L r 


2.5 


0.8 


0.8 


1.3 


2.6 


IRS spectrum 


Weedman et al. 2006 (15) 




Weedman et al. 2006 (9) 


Weedman et al. 2006 (24) 


Weedman et al. 2006 (5 


IRS features 


No Si Abs. 




Possible Si Abs. 


Possible Si Abs. 


featureless 



a Abbreviated from the full name of SST24 J142827.1+354127 etc. 

b Note that SST24 J14301 1.3+343439 is not strictly an optically faint ULIRG by our definition of J?- [24] > 14. 

c The quoted positions of our targets are measured directly from the optical images. In cases where the source is optically invisible, we used the MIPS position corrected to be on the 
NDWFS frame. 

d In cases in which no catalog magnitude exists, we quote the magnitude obtained within a 4 "aperture in brackets. 

e For both the IRS and optical redshifts, - denotes that the object was not observed and ? denotes that a redshift could not be determined. The IRS redshifts are accurate to within 0.1. 
f The optical redshifts were obtained on Keck LRIS and DEIMOS (Desai et al. in preparation). 

6 The continuum amplitude is given in units of 1E-18 erg s _1 cirT 2 A -1 and refers to the fitted value in around the Ha/[NII] or H/37[OIII] lines, 
k All flux densities are given in units of 1E-18 erg s _1 cm" 2 . 

' All values refer to [Nil] A6583. The [Nil] A6548 a is fixed to be the same as that of [Nil] A6583 and the [Nil] A6548 flux is fixed to be 1/3 of [Nil] A6583. 
J All values refer to [OIII] A5007. The [OIII] A4959 <r is fixed to be the same as that of [OIII] A5007 and the [OIII] A4959 flux is fixed to be 1/3 of [OIII] A5007. 
k The luminosities are from the measured flux densities and do not include any corrections for the substantial correction due to slit losses (see Section 3.4) or dust extinction (see 
Section 5.3). 

' Sources are classified as AGNl if they exhibit broad (>1000 km s~') Balmer lines and as AGN2 if they exhibit narrow Balmer lines and [NII]/Ha and [OIII]/H/3 ratios characteristic 
of AGN activity. 

m The intrinsic infrared luminosity, Lir(8- 1000 <jm), is estimated from L24 using a bolometric correction of 7 (typical for AGN-dominated ULIRGs; Xu et al. 2001) 
" The Gaussian ler width of the emission line (=2.263*FWHM) 
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from optical spectroscopy, was targeted in the Zf-band to de- 
termine line strengths and widths of the [OIII] AA5007, 4959 
doublet and H/3 lines. Exposure times are listed in Table 1; 
they were typically wl hour per filter. 

3.2. Gemini /NIRI 

Six sources were observed with NIRI on the Gemini tele- 
scope on 7th - 8th June 2006 UT, and 29th July 2006 UT in 
queue mode. We used the NIRI f / 6 with a grism with the 4- 
pix (« 0"4 wide) "blue" slit. All 6 sources that were observed 
were detected in at least one of the J (G0202), H (G0203), and 
K (G0204) bands. Four sources exhibited continuum emis- 
sion but no obvious emission lines or breaks. Two sources 
exhibited strong emission lines and form part of the sample 
that is presented in this paper. SST24 J143027. 1+344007 
was targeted in the J (G0202) and H (G0203) bands; SST24 
J 142800.7+350455 was targeted in the H (G0203) and K 
(G0204) bands (see Table 1). These grisms correspond to res- 
olutions of R«650, 940, and 780 for the J, H, and K filters 
respectively. For sky subtraction, the observations were split 
into 4 minute exposures, and the object was nodded along the 
slit by 10". Typical exposure times were «1 hour per filter 
(see Table 1). 

3.3. Data reduction 

The OH telluric sky lines are many orders of magnitude 
brighter than the signal from our science targets and must be 
accurately subtracted. We began by removing the first-order 
background spectrum by subtracting the closest observation in 
time with a different dither position. Both the NIRSPEC and 
NIRI low-resolution spectra are tilted relative to the array and 
are not aligned with the detector columns or rows (McLean 
et al. 1998). We corrected for this by rectifying the spectra. 
This requires some interpolation between adjacent pixels but 
makes the fitting of the residual background and extraction 
of the spectra much simpler. Any residual background sky 
lines (resulting from the time-variability of telluric sky lines) 
were then removed by fitting an 15th order polynomial along 
the lines (masking out the position of the target and any other 
objects in each frame). The 2-D images were then shifted to 
move the object spectrum to the same column, and averaged 
to obtain a final 2-D image. Sigma clipping was used to reject 
cosmic rays and a rectified bad pixel mask (obtained from the 
median of multiple dark exposures) was used to reject bad 
pixels during stacking. The 1-D spectrum was then extracted 
from the stacked 2-D image and wavelength calibrated using 
the night sky lines and Argon arcs. All our reductions used 
standard tasks in the IRAF LONGSLIT package. 

3.4. Flux calibration 

Throughout our observing run, we observed telluric stan- 
dard stars at a similar airmass to our target observations in 
order to correct for telluric absorption and provide relative 
flux calibrations of the spectra. We used F0V stars as a com- 
promise between avoiding strong hydrogen lines which are 
stronger in hotter stars and the lines of other atomic species 
which become noticeable in cooler stars. For flux calibration, 
we used the following procedure. We obtained a model spec- 
trum of an F0V star scaled to the magnitude appropriate for 
the standard star (using the IRAF task PLSPEC). We then di- 
vided this by the observed spectrum of the standard star to 
create a calibration spectrum, which we multiplied into the 
source spectrum. 



We note that the flux calibration may have large uncer- 
tainties due to slit losses. Although the seeing was gener- 
ally good («0."5-0"7), the slit was only 0"76 and 0."4 wide 
for NIRSPEC and NIRI respectively and the targets were too 
faint to position on the slit precisely. Any offset in the tar- 
get with respect to the slit will have resulted in an underes- 
timate of the flux density of the source. For 4 sources, we 
have TT-band detections from the NDWFS or FLAMEX sur- 
veys or from observations by Keck NIRC. In the 2 cases with 
a significant TT-band continuum detection, we compare the K- 
band magnitudes with that obtained from the spectra to check 
our calibrations. For SST24 J142644.3+333051 and SST24 
J143424.4+334543, the flux is 0.78 and 0.64 of that expected 
from their TT-band magnitudes. For consistency with the ob- 
servations with no TT-band magnitudes, we chose not to cor- 
rect the flux. However, we note that the absolute flux mea- 
surements are likely to be underestimated by «25-50%. 

3.5. Emission line fitting 

We determined line widths and strengths by fitting Gaussian 
models to the 1-D spectra. These lines are fit simultaneously 
using MPFITFUN: a robust non-linear least squares curve fit- 
ting routine. Each line is approximated by a Gaussian profile 
with the width, height and position allowed to vary. The con- 
tinuum flux is assumed to be constant as a function of wave- 
length and the height is simultaneously fit with the height, 
width, and positions of the Gaussian profiles. In order to min- 
imize the number of fitted parameters and help fit the nois- 
ier lines, we used the following method. We tied the central 
wavelengths of the individual lines to have the same redshift. 
We only allowed the fitted redshift to differ from the redshift 
measured by eye by 0.03 and constrained the fits to only re- 
turn positive emission line strengths. We also required all the 
forbidden lines to have the same widths. We fixed the [Nil] 
A6583 / A6548 ratio to 3 (as derived by Storey & Zeippen 
2000) for the transition probability ratio since these lines are 
often blended with Ha and typically have low signal-to-noise 
ratios. We also fixed the [OIII] A5007 / [OIII] A4959 ratio to 
its theoretical value. We chose not to fit the [SII] AA 6717, 
673 1 lines since, if they are present at all, they have very low 
signal-to-noise ratios. In cases for which we had measure- 
ments of strong [OIII] lines and poor constraints in the [Nil] 
lines, we fixed the width of the [Nil] lines to that of the [OIII] 
lines. The line fitting was weighted by the inverse variance 
spectrum. This gives less weight to wavelengths which are 
affected by strong (and often hard to perfectly subtract) sky 
lines. Figure 1 shows an example of how the fitting process 
works. 

4. RESULTS 

4.1. The near-IR spectra 

The final flux-calibrated near-IR spectra are presented in 
Figure 2. The fitted line profiles are overplotted. In all 9 spec- 
tra covering the appropriate rest wavelengths, Ha is promi- 
nent and has a highly significant detection over the noise. It is 
broad (> 1900 km s _1 ) in 6 cases (we assume that the broad- 
line component in SST24 J142827. 1+354127 dominates the 
24^m luminosity; see Section 4.2). [Nil] AA6548,6583 and 
[SII] AA67 17,6731 are also present in many cases, although 
they are not always significant detections. In the 5 spectra 
covering the rest wavelengths of H/3 and [OIII] AA4959,5007, 
[OIII] A5007 has a significant detection in all but one case, 
[OIII] A4959 has a significant detection in two cases, and 
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FIG. 1 .— Example of line fitting for SST24 J143028.5+343221. The blue 
histogram shows the un-smoothed data. The model fit is shown by the dashed 
red line. We also show the separate contributions to the model fit from the 
continuum and individual Gaussian line profiles (black dotted lines). Wave- 
length regions particularly affected by sky emission lines are shaded in yel- 
low. These regions are given less weight in the fitting process. 

H/3 has only a marginal detection in one case. We present 
the fitted widths and fluxes of the observed emission lines 
in Tables 1. All quoted line widths and fluxes have been 
corrected to take into account the instrumental profile of 
FWHM,„ J( wlOA and «14A for the H- and /T-band filters 
respectively on NIRSPEC and FWHM,„ S , «24A, «20A, and 
«28A for the /— , H-, and /f-band filters respectively on 
NIRI, as measured from the width of the night sky lines. We 
discuss the individual sources (on order of increasing redshift) 
in the following section. 

4.2. Notes on individual sources 

SST24 J142827.1 +354127.- -This source is highly unusual in 
its powerful 24 /im emission (f24=10 mJy) and faint (but very 
red) optical magnitude (Bw-K=ll). The source is discussed 
in detail in Desai et al. (2006), where both its IRS and optical 
spectra are presented. The optical spectrum exhibits a single 
emission line which is identified as [Oil] providing a redshift 
of z =1.292. Given this redshift, SST24 J142827. 1+354127 
has one of the largest rest-frame 5 /im luminosities known 
(z/L„ = 5.5 xlO I2 L©), yet its faint optical luminosities and 
lack of X-ray emission down to the limits of the XBotites 
survey imply extreme gas and dust obscuration. The IRS 
spectrum is typical of that of AGN-dominated sources: it 
exhibits a strong and silicate absorption feature but no PAH 
features. Given its known redshift, we observed SST24 
J142827. 1+354127 in the //-band to target the Ha line. Our 
near-IR spectrum shows significant continuum and inter- 
estingly, reveals two redshift components which we denote 
M142827a and M142827b. The redshift of M142827a agrees 
with that of the optical spectrum. Unfortunately, a powerful 
night sky line happens to fall at the same wavelength as the 
observed wavelengths of the diagnostically important lines, 
Ha and [Nil] A6583, making estimates of the Ha flux and 
width very uncertain. Given the effect of the sky line, our 
determination of FWHMh q ~230 km s -1 is highly uncertain. 
However, it is unlikely to be broad enough to constitute a 
broad AGN emission line. The [NII]/Ha ratio is poorly 
constrained and so we classify this component as dominated 
by either a starburst or type 2 AGN. M142827b is newly 
discovered in our NIRSPEC observations. It is revealed by 
a broad line at an observed wavelength of 15,689 A. If this 



is Ha, then the source is at z =1.383. Given the width of 
the line, the source must be dominated by a type 1 AGN. 
Although there is no evidence for a second source in the 
optical spectrum, the IRS redshift estimated from the silicate 
absorption feature falls between the two NIRSPEC redshifts 
as one would expect for a composite source in which both 
components were optically faint. Given the large velocity 
separation between the two redshift components (Av~l 1,500 
km s" 1 ), we conclude that the mid-IR flux is coming from two 
separate sources which happen to be coincident along our line 
of sight. If SST24 J142827.1+354127 is a composite source, 
this would certainly explain its unusually high infrared flux. 
Because the Ha luminosity of the broad line is «10 times 
larger than that of the narrow Ha line of the lower redshift 
source, we assume that M142827b is the dominant source 
of the infrared flux, and hence the source which satisfies our 
R-[2A] > 14 criterion. 

SST24 J 143027. 1+344007.— This source was observed 
with Gemini NIRI because of its bright /f-band magnitude 
(/sT=17.89). It exhibits broad Ha emission (FWHM„ Q =1940 
km s" 1 ) in the //-band confirming that it is dominated by an 
AGN, and providing a redshift of z=1.370. Our observations 
in the /-band show a marginal detection of [OIII]A5007 and 
no H/3 detection. The source exhibits significant continuum 
emission and so cannot have simply been missed off the slit. 
The weakness of [OIII] emission is unusual in comparison to 
the other sources in our sample. Using the minimum Ha to 
H/3 ratio, we estimate at least 4.6 magnitudes of extinction at 
the wavelength of Ha. The unattenuated Ha luminosity must 
be at least 3 x 10 44 ergs s" 1 . 

SST24 J143312.7+34201L— This source was observed in 
the ZT-band by Keck NIRSPEC because its redshift (z = 2.2; 
as estimated from the silicate absorption feature) put the 
observed wavelength of Ha in the /if -band. Apart from a 
moderate silicate absorption feature, the IRS spectrum is 
featureless. It is best fit by an AGN-dominated template. Our 
NIRSPEC observations show that the source has a slightly 
lower redshift of z =2.114, putting Ha at the very edge of 
our wavelength coverage. We do however, appear to have 
observed the peak of the line, so our estimates of the line 
strength and width should not be too severely affected. The 
Ha line is very broad (FWHM H «=3150±300 km s" 1 ) with 
a large equivalent width confirming that the source is domi- 
nated by an AGN. It has the second brightest Ha luminosity 
of our sample. The presence of [Nil] is not significant given 
the large uncertainties introduced by a coincident night sky 
line. The source appears to have [SII] line emission, although 
at fairly low significance. 

SST24 J14301 1.3+343439.— This source was chosen 
because its redshift from Keck optical spectroscopy (z=2.12) 
put the observed wavelength of Ha in the Zf-band. It was 
observed in both the K- and //-bands. With R- [24] = 12.5, 
this source does not strictly obey our R- [24] > 14 selection 
criteria and it has not been observed with IRS. The near-IR 
spectrum confirms the optical redshift (z=2.114). It exhibits 
relatively narrow Ha (FWHM=470 km s" 1 ) and its combined 
[NII]/Ha ratio and [OIII]/H/3 limits suggests that it is a type-2 
AGN. From the minimum Ha to H(3 flux ratio, we estimate 
at least 2.4 magnitudes of extinction at the wavelength of Ha. 
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FIG. 2.— Keck NIRSPEC and Gemini NIRI near-IR spectra of the 10 optically faint ULIRGs. The spectra have been smoothed by a 14A, 22A, 10A, 15A, 
and 35A boxcar filter in the NIRSPEC H- and K— band filters and NIRI J, H, and K filters respectively (solid line). The expected positions of typical emission 
lines are labeled. For each object, the thin solid line (offset from zero by up to 2 X 10~' 8 erg s cirT 2 A~' for clarity) shows the corresponding 1 a error spectrum 
which is dominated by the OH telluric emission lines. Regions that are particularly affected by the sky emission lines are shaded in yellow. The two-dimensional 
spectra are also shown with the same wavelength range as that of their 1-D spectra. The fitted line profiles are shown by the dashed line. In the 4 cases in which 
the H/3 line is not detected, the 3<r minimum line profile (with the width fixed to that of Ha) are shown by the dotted line. 
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SST24 J143028.5+34322L— This source was originally 
chosen for near-IR spectroscopy because optical follow-up 
on Keck LRIS found a redshift of z = 2.178 (from Lya). It 
was then subsequently observed with Spitzer IRS which con- 
firmed this redshift. The source exhibits relatively narrow Ha 
(FWHM=542 km s" 1 ) and its large [NII]/Ha ratio suggests 
it is a type-2 AGN. The IRS spectrum shows a deep silicate 
absorption feature and possible PAH emission features (see 
Hidgon et al. in preparation for more details). 

SST24 J 142842.9+342409.— This source was chosen 
because it exhibited a featureless power-law spectrum in the 
mid-IR. Follow-up optical spectroscopy also failed to yield 
a redshift. More sophisticated fitting of the IRS spectrum in 
Weedman et al. (2006) estimated a possible redshift of z =1.1. 
We observed it in the /f-band and detected a strong Ha line, 
identifying the source as AGN-dominated and putting it at a 
redshift of z=2.180. We then followed up in the //-band to 
target the [OIII] and H/3 lines. The Ha line is broad («3000 
km s" 1 ) with a large equivalent width and the observed Ha 
luminosity of w 2 x 10 ergs s" 1 is the largest in the sample. 
The [OIII] line is unusually broad (FWHM [o///] =1600 
km s" 1 ), perhaps due to shock-induced superwind activity 
(cf. Takata et al. 2006). The //-band spectrum shows no 
significant H/3 flux suggesting large columns of attenuating 
dust (see Section 5.3). Given the new redshift, we re-examine 
the IRS spectrum but find no evidence for any silicate 
absorption in the source. This is somewhat surprising given 
the large columns of dust which must exist in this source. See 
Section 5.4 for further discussion of this source. 

SST24 J "142800.7 ' +350455.- -This source was observed 
with Gemini NIRI because of its bright ZT-band magnitude 
(Zf=17.78). It has not been observed with IRS. The strong 
[OIII] and Ha lines detected in the H— and Zf-bands respec- 
tively determine that z=2.223. Although the best obtained fit 
gives FWHMtf Q = 3833 km s" 1 , other (slightly worse) fits 
indicate a narrower line profile. Nevertheless, all fits suggest 
FWHM# a >2000 km s" 1 and the source is classified as a type 
I AGN. Because the Ha flux is relatively low and FWHM// Q 
is relatively high, we can only put weak constraints on the 
extinction. 

SST24 J143424.4+334543.— This source was chosen 
because, although the IRS spectrum showed a featureless 
power-law spectrum, optical follow-up spectroscopy found a 
redshift of z = 2.26. The NIRSPEC spectrum confirms the op- 
tical redshift and shows a very large [NII]/Ha ratio. Although 
difficult to measure given the powerful [Nil] emission, the 
width of the Ha line is relatively narrow and we classify it as 
a type 2 AGN. The IRS spectrum shows a downward trend 
in the flux at «31/im where one would expect the silicate 
absorption feature to lie. It may be that the feature is too near 
the edge of the spectrum to yield a convincing IRS redshift. 
The IRS spectrum shows possible PAH emission features al- 
though it is fairly noisy. It is interesting that the only sources 
with IRS spectra showing possible PAH emission features 
are the 2 sources which are optically classified as type 2 AGN. 

SST24 J 142939.1 +353558.— -This source was observed 
with IRS and Keck LRIS but no redshift was determined. 
We obtain a redshift of z=2.498. The Ha line is broad 
and we classify the source as a type 1 AGN. The source is 



also detected in the Chandra XBootes survey (Kenter et al. 
2005; Murray et al. 2005; Brand et al. 2006a) with an X-ray 
luminosity of Lo.5-7^v=7.7x 10 42 ergs s" 1 . The IRS spectrum 
did not yield a redshift because the silicate absorption feature 
falls at the very edge of the spectrum. 

SST24 J142644.3+333051.— We targeted this source 
because of its featureless IRS spectrum. We observed it in 
the ZT-band to target the H/3 and [OIII] lines, given its optical 
redshift of z = 3.355 obtained on Keck LRIS (see Desai et 
al. in preparation). We obtain a slightly lower redshift of 
z= 3.312 from the [OIII] lines. Although only marginally 
detected, the H/3 line is broad (FWHM =2066±394 km 
s" 1 ) and we classify this source as a type 1 AGN. SST24 
J142644. 3+333051 is also a powerful X-ray source with 
Lo.5-7fev=9x 10 44 ergs s" 1 , confirming that it must be a 
powerful AGN. The silicate absorption feature falls beyond 
the observed wavelength range of the IRS spectrum. 

5. THE CHARACTERISTICS OF OPTICALLY FAINT ULIRGS 

In this section, we use line diagnostics to characterize the 
sources into broad-line type 1 AGN, narrow-line type 2 AGN 
or starburst-dominated sources. We derive physical properties 
and in cases with both Ha line measurements and upper lim- 
its for the H/3 line flux, we estimate lower limits to the dust 
extinction. We also discuss the mid-IR spectra. 

5.1. Line diagnostics 

Six sources show broad (>1900 km s" 1 ) Ha lines which 
are characteristic of AGN-dominated sources. For the one 
source observed at shorter rest-frame wavelengths (SST24 
J142644.3+333051), the H/3 line is only marginally detected 
but also appears to be broad. We also classify this source 
as a type 1 AGN. A common diagnostic to distinguish be- 
tween AGN- and starburst-dominated sources uses the line 
ratios of [NII]/Ha and [OIII]/H/3 as indicators of the hard- 
ness of the EUV radiation field in the narrow-line regions 
(Baldwin et al. 1981; Osterbrock & Pogge 1985; Veilleux 
& Osterbrock 1987; Kauffmann et al. 2003; Kewley et al. 
2006). Unfortunately, these diagnostics are unreliable for 
sources with broad Balmer lines since the broad-line compo- 
nent, which likely arises in the broad-line region (BLR), typ- 
ically overwhelms any narrow-line component. Three of our 
sources show narrow-line Ha emission and have significant 
[Nil] emission, enabling reliable measurements of the ratio. 
The [NII]/Ha ratios are 0.33, 0.48, and 4.6. AGN-dominated 
sources typically have [NII]/Ha >0.7 (e.g., Swinbank et al. 
2004) suggesting 1 of our sources is a type 2 AGN. The source 
with [NII]/Ha=0.33 has a high [OIII]/H^ ratio. The [NII]/Ha 
vs. [OHTJ/H/3 classification schemes of Kewley et al. (2001) 
and Kauffmann et al. (2003) suggest that this source is also a 
type 2 AGN. Our limited statistics therefore suggest that both 
the broad and narrow-line sources in our sample are AGN- 
dominated. 

5.2. Physical properties 

The main derived physical properties of the sources are 
shown in Table 1. The sources have Ha luminosities in the 
range 1 x 10 42 -2.4 x 10 43 ergs s" 1 . For sources classified as 
type 1 AGN, the SMBH mass can be estimated from both 
the FWHM and luminosity of the Ha line using the prescrip- 
tion of Greene & Ho (2005). We estimate values of 0.2-1.3 
x 10 8 M Q which are surprisingly low for such bolometrically 
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luminous sources. However, as we discuss in Section 5.3, 
these values are likely to be severely underestimated due to 
attenuation by dust. To determine the infrared luminosities 
from the observed 24 pm luminosities, we assume a bolo- 
metric correction of 7.0 (typical of AGN-dominated sources; 
Xu et al. 2001). All our sources have L 1R > 1O 12 L and are 
therefore classified as ULIRGs. Two sources are detected 
in the Chandra XBootes survey (SST24 J142939. 1+353558 
and SST24 J142644.3+333051 with X-ray luminosities of 
Lo.5-7fcv=7.7x 10 42 ergs s" 1 and 9xl0 44 ergs s" 1 respectively. 
For SST24 J142939. 1+353558, assuming a canonical accre- 
tion efficiency of e=0. 1 and that the correction from the X-ray 
luminosity to the bolometric luminosity is a factor of 35, the 
accretion rate onto the SMBH is « 5 x 10 7 M Q Gyr _1 (Barger 
et al. 2001). Assuming the Ha line is attenuated by a factor 
of 10 in flux, its black hole mass is « 10 8 M Q . If the source 
has been accreting at a similar rate throughout its lifetime, it 
would have taken it ~ 2 Gyrs to build up to this mass. Given 
this large value, we suggest that it is likely that the X-ray lu- 
minosity is also attenuated by large columns of gas associated 
with the dust. This is also supported by the fact that only 2/10 
of the sources are detected in the XBootes survey even though 
they are classified as powerful AGN. 

5.3. Dust extinction 

Observations of the Balmer decrement can provide an esti- 
mate of the extinction due to dust in these sources and hence 
the correction one must apply to obtain intrinsic line lumi- 
nosities. Gaskell & Ferland (1984) show theoretically that 
for AGN, Ha/H/3 «3.1. This is slightly higher than the ex- 
pected value of Ha/H/3=2.86 for case B recombination at 
T=10,000Kandn e «10 4 cm" 3 (Osterbrock 1989). We assume 
Ha/H/3=3.1 for the intrinsic (unreddened) ratio. In three cases 
(SST24 J143027.1+344007, SST24 J14301 1.3+343439, and 
SST24 J 142842.9+342409), we have significant detections of 
the Ha line flux and useful upper limits for the H/3 line flux. 
We measure 3a upper limits on the H/3 flux from the noise 
spectrum assuming the same FWHM as measured for the Ha 
line (see Figure 2). The upper limits are given in Table 1 
and correspond to Ha/H/3 > 22.5, 8.5, and 15.7 for SST24 
J143027.1+344007, SST24 J14301 1.3+343439, and SST24 
J 142842.9+342409 respectively. This demonstrates that both 
the broad- and narrow-line sources suffer from significant op- 
tical extinction. Assuming the standard extinction curve of 
Osterbrock (1989), we obtain E(B-V)> 1.9, 1.0, and 1.6, cor- 
responding to extinctions of A(H Q )>4.6, 2.4, and 3.8 magni- 
tudes. This corresponds to a correction factor of at least , 69, 
9, and 33 times the Ha luminosities quoted in Table 1. The 
extinction limits are far larger than that inferred for composite 
spectra of reddened quasars from SDSS (Richards et al. 2003) 
and larger than even the largest extinctions inferred for UV- 
selected galaxies (A(H Q )=1.7; Erb et al. 2006) and distant red 
galaxies (A(H Q )=1.7; van Dokkum et al. 2004) at similar red- 
shifts. If the true extinctions are close to the measured limits, 
although they sample sizes are small, they may be compara- 
ble to sub-millimeter galaxies (Takata et al. 2006), local dusty 
starbursts (A(H Q )=2.9; Poggianti & Wu 2000, and some of the 
largest optical extinctions of the most extreme local ULIRGs 
by Veilleux et al. (1995). Assuming that all the sources suffer 
a similar extinction, the line widths and corrected Ha lumi- 
nosities of > 10 43 - 10 44 ergs s" 1 suggest SMBH masses of 
> 10 8 -10 9 M Q . 

In Figure 3, we show the infrared luminosity versus the 
measured Ha luminosity for our sample. Our sources have 



much smaller Ha luminosities relative to their infrared lumi- 
nosities when compared to local IRAS galaxies (Kewley et al. 
2002). We compare our results to the sub-mm selected sam- 
ple of Swinbank et al. (2004). Their sample provides a good 
comparison sample given the similar redshift range and in- 
frared luminosities. Their observations were also taken with 
the KECK / NIRSPEC 0.76" slit, so slit losses should be simi- 
lar. A Kolmogorov-Smirnov test (Press et al. 1992) shows that 
there is only 5% chance that the Ha luminosities are drawn 
from the same distribution. This suggests that optically faint, 
infrared bright galaxies tend to have higher optical extinctions 
than sub-mm selected galaxies. 
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FIG. 3. — The observed Ha luminosity versus the infrared luminosity for 
the 9 sources with Ha line measurements (large filled circles). For compar- 
ison, we plot a sample of local IRAS galaxies (Kewley et al. 2002; small 
empty circles), and sub-mm selected galaxies (Swinbank et al. 2004; open 
squares). The dotted line shows where the infrared and Ha SFR indicators 
would agree. The Ha luminosities are smaller than we would expect given 
the large infrared luminosities, suggesting large optical extinctions. 



5.4. Comparision to mid-IR spectra 

In Houck et al. (2005) and Weedman et al. (2006), we pre- 
sented mid-IR spectroscopy of R- [24] > 14 sources using 
Spitzer IRS. From fitting of the IRS spectra with templates, 
Houck et al. (2005) conclude that AGN probably dominate 
the mid-IR emission in the majority of > 0.8 mJy op- 
tically faint ULIRGs. Our results are consistent with this. 
Weedman et al. (2006) presented the spectra for which we 
had previously been unable to determine a confident redshift 
in Houck et al. (2005). 7/10 of our NIRSPEC sample were 
observed with Spitzer IRS as part of our larger IRS program. 
Three sources show strong silicate absorption features in their 
IRS spectra and have IRS redshift determinations that are 
broadly consistent with their NIRSPEC redshifts. The four 
sources with no clear spectral features in their IRS mid-IR 
spectra have redshifts derived from their near-IR spectra of 
z=2.180, 2.263, 2.498, and 3.312 (see Table 1). The 9.7 ^m 
silicate absorption feature is redshifted beyond the IRS spec- 
tral window for redshifts z > 2.5, thus naturally explaining the 
featureless IRS spectra for the two highest redshift sources. 
The IRS spectra of SST24 J143424.4+334543 at z=2.263 and 
SST24 J142939. 1+353558 at z=2.498 do show downturns in 
the flux at the long wavelength edge of the IRS window, sug- 
gesting the onset of the silicate absorption. However, SST24 
J 142842.9+342409 at z=2.180 is at low enough redshift to 
have shown any silicate absorption feature; its absence sug- 
gests a perhaps unusually weak feature, and we discuss this 
further in Section 6.3. 
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We also observed two further IRS sources (SST24 
J142611.3+351217andSST24J142850.9+353146: source 20 
and 23 in Weedman et al. (2006)) with NIRSPEC but detected 
no significant continuum or line emission. The Zf-band mag- 
nitudes are such that we may not be able to detect continuum 
emission. Assuming that the source did not fall out of the 
slit, it is likely that these sources are at redshifts for which 
none of the powerful lines fell into the Zf-band spectral cover- 
age. SST24 J14261 1+351217 is predicted by Weedman et al. 
(2006) to be at z w 1 .6 which is too low for us to have detected 
Ha in the /sT-band. SST24 J 142850.9+353 146 has a feature- 
less mid-IR spectrum. If it is in the redshift range 2.7< z <3.1 
or at z >3.7, none of the typically strong lines would have 
been detectable. 

6. DISCUSSION 

6.1. The number density of obscured AGN 

Our near-IR spectra suggest that most f24 >0.8 mJy opti- 
cally faint luminous infrared galaxies harbor AGN. This is an 
important confirmation of our Spitzer IRS results and shows 
that there is a significant population of AGN missing from tra- 
ditional optical surveys. Even giving conservative estimates 
of their bolometric luminosities, these sources are comparable 
in luminosity to bright optically selected quasars. Brown et al. 
(2006) found that the number density of 24 /im selected, op- 
tically unobscured AGN was comparable to the number den- 
sity of optically selected AGN. They find «140 f24 >1 mJy 
quasars with z > 1.3 in the entire Bootes field. There are 
«340 f24 >1 mJy sources with R— [24] >14. If we assume 
that these sources are all at z > 1.3 and harbor luminous AGN 
(as is suggested by our results and their large bolometric lu- 
minosities), then they constitute an important population of 
AGN with space densities of at least twice that of the opti- 
cally luminous type I AGN. Although this is only a crude es- 
timate, it agrees broadly with previous estimates of the num- 
ber of obscured to un-obscured AGN (e.g., Ueda et al. 2003; 
Martmez-Sansigre et al. 2005). 

6.2. Obscuration by dust 

7/10 of the sample presented in this paper are broad-line 
(optical type I) AGN. This fraction agrees broadly with re- 
ceding torus models which predict that the broad-line region 
is surrounded by an optically thick molecular torus (e.g., An- 
tonucci 1993) and that the opening angle of the torus increases 
with luminosity because more luminous AGN can sublimate 
the dust out to larger distances (Lawrence 1991; Arshakian 
2005; Simpson 2005). Simpson (2005) predict the fraction of 
type I AGN as a function of [OIII] luminosity. Our typical 
intrinsic [OIII] luminosities of w 10 43 — 10 45 ergs s" 1 imply 
type I fractions of 0.6-1 depending on the adopted model. On 
first look, the fraction of broad line sources is consistent with 
what we would expect for such luminous AGN. However, this 
population was selected to have faint optical (i.e., rest-frame 
UV) magnitudes. Given that all of the sources have R>22, 
we might expect our sample to be biased to heavily extincted 
quasars, and therefore only visible as type II AGN. The visi- 
bility of the broad lines suggests that the lines of sight to the 
broad line region are not completely obscured and hence that 
it is unlikely that these quasars are being viewed through the 
mid-plane of a dusty, completely optically thick torus. Our 
results are also consistent with models in which the torus is 
identified with the dusty, optically thick region of the wind 
coming off the central accretion disk (e.g., Elvis 2006; Elitzur 



& Shlosman 2006). If the obscuring clouds are clumpy, this 
may allow one to view the central broad-line regions in some 
cases. 

The Ha / H/3 Balmer decrements are large for both narrow- 
line and broad-line AGN in our sample. This, along with 
the optical faintness, suggests that the AGN continuum emis- 
sion, broad-line region, and narrow-line regions are heav- 
ily extincted. Although the extinction of the AGN contin- 
uum and broad-line regions only requires extinction on small 
scales, extinction of the narrow-line region requires extinc- 
tion on much larger scales. We estimate SMBH masses of 
~ 10 8_9 Mq. Given the large infrared luminosity of the AGN, 
the central SMBH must also be undergoing very rapid growth. 
If the Magorrian relation holds, we would expect the host 
galaxies to have large bulges with significant accompanying 
starbursts which should be visible in the rest-frame UV. We 
see much fainter optical emission than we would expect, sug- 
gesting obscuration on very large (i.e., kpc) scales. The fact 
that we see hints of extended but very faint optical emission in 
many of our sources supports this view (Dey et al. in prepara- 
tion). 

Martmez-Sansigre et al. (2006) propose two types of ob- 
scured AGN. The first are the traditional type II AGN whose 
optical light is obscured by a dusty torus ("torus-obscured" 
AGN). The second are AGN whose optical light is obscured 
on much larger scales ("host-obscured" AGN). They suggest 
that although receding torus models predict type II to type I 
ratios of ~ 1 : 1 , larger ratios may be needed to explain the hard 
X-ray background (~ 3 -4 : 1 ; Gilli et al. 2001). The discrep- 
ancy could be solved by the additional population of "host- 
obscured" AGN which will be counted as obscured AGN in 
the latter case but exhibit both type-I and type-II optical spec- 
tra. As we discuss above, the selection of optically faint lu- 
minous infrared galaxies will naturally select "host-obscured" 
sources. Albeit with only a small sample, our near-IR spectra 
show that these sources indeed have the fraction of broad- to 
narrow-line sources that one would expect for such luminous 
sources, despite being too obscured to appear in optically se- 
lected samples. 

6.3. Silicate absorption and optical extinction 

There is no obvious correlation between the silicate absorp- 
tion depths in the IRS spectra and optical rest-frame proper- 
ties. SST24 J1433 12.7+34201 1 exhibits deep silicate absorp- 
tion in its IRS spectrum whereas SST24 J142842.9+342409 
must have only shallow silicate absorption yet both have 
very similar rest-frame optical spectra (with prominent broad 
Ha emission lines). SST24 J1433 12.7+34201 1 and SST24 
J143028. 5+343221 both have deep silicate absorption fea- 
tures yet they exhibit broad and narrow Ha emission lines 
respectively. 

Are the silicate depths consistent with the observation of 
broad emission lines? SST24 J1433 12.7+34201 1 has broad 
emission lines and has relatively deep silicate absorption 
(A(9.7/im)^1.2). Assuming the dust models of Li & Draine 
(2001) which predict a extinction ratios of A(H Q )/A(9.7^m) 
«10, this corresponds to A(H a )«12 which should be too 
large for us to have detected broad Ha for any plausible in- 
trinsic Ha luminosities (L(Ha) < 10 46 ergs s" 1 ). However, this 
apparent discrepancy also exists for some local ULIRGs. Hao 
et al. (2006) investigate the distribution of silicate absorption 
strengths for a wide variety of local AGN and ULIRGs. For 
type I AGN, the deepest silicate absorption strength is found 
for Mkn 231 (Weedman et al. 2005) with A(9.7/zm)=0.75 
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which would correspond to A(H a )«7.5, again implausibly 
high. Arp 220 has deep silicate absorption and evidence for 
a heavily embedded nuclear source (e.g., Spoon et al. 2004), 
yet optical emission lines are visible and lead to a liner classi- 
fication. We suggest that in these cases, the dust distribution 
may have a complicated geometry. Possible scenarios include 
a less than 100% covering factor due to a clumpy dust distri- 
bution or scattering of the optical light. If the sources have op- 
tically thick nuclei, the situation may be further complicated 
by the fact that we may be seeing down to different optical 
depths at different wavelengths. 

SST24 J 142842.9+342409 is an interesting case given its 
lack of any silicate absorption feature. The silicate absorp- 
tion feature should be observable at an observed wavelength 
of «31 fim but no such feature is observed (see Weedman 
et al. (2006) figure 1; source 15). This is particularly interest- 
ing given that Ha/H/3 > 15.7, corresponding to E(B-V)>1.6 
and suggesting large columns of attenuating dust. Could the 
geometry of this source result in a weaker silicate absorp- 
tion feature? The presence of broad Ha suggests that this 
source is a type 1 AGN in which the torus is viewed roughly 
face-on. Type 1 AGN are expected to show silicate emis- 
sion from the hot illuminated surface of the inner torus (Pier 
& Krolik 1992). For more than a decade, silicate in emis- 
sion in type 1 AGN remained elusive. However, they have 
recently been found in powerful unobscured AGN (Sieben- 
morgen et al. 2005; Hao et al. 2005) and in the ULIRG FSC 
10214+4724 (Teplitz et al. 2006). Siebenmorgen et al. (2005) 
suggest that this may be a luminosity effect, with only the 
most luminous AGN showing silicate in emission. If this 
is the case, then we might expect to see silicate emission in 
our sample of very luminous type 1 optically faint ULIRGS. 
SST24 J 142842.9+342409 has the second largest infrared lu- 
minosity of sources at low enough redshifts to observe the 
silicate absorption feature. Perhaps the combination of sili- 
cate emission from the torus and silicate absorption further 
out, results in only a weak silicate feature in this case. 

7. CONCLUSIONS 

We have obtained near-infrared spectroscopy for a small 
sample of 24 fim sources with extremely red optical-to-mid- 
infrared colors selected from the Spitzer MIPS survey of the 
NDWFS Bootes field. We conclude that: 

(1) All of the sources in the sample lie at high redshift 
(1.3 < z < 3.4) and are therefore very luminous infrared 
sources. 

(2) All the sources show either broad-lines or have narrow- 
line ratios suggesting that they contain powerful AGN. 
This is an important confirmation of our Spitzer IRS 
results and shows that there is a significant population 
of AGN missing from traditional optical surveys. If all 
R- [24] > 14 sources with f 24 > 1 mJy harbor AGN, 
then these sources may be twice as common as more 
optically luminous type I AGN at similar redshifts. 

(3) 7/10 of the sample have broad (> 1900 km s" 1 ) Ha or 
H/3 emission lines. We classify these sources as type I 
AGN. The remaining 3 sources have Balmer line widths 
> 450 km s" 1 and line diagnostics predict them to be 
type II AGN. 

(4) The fraction of type I to type II AGN is consistent with 
receding torus models given their large intrinsic [OIII] 



luminosities. This is somewhat surprising: given their 
high infrared luminosities and faint optical magnitudes, 
we might expect these sources to all be heavily ex- 
tincted quasars (i.e., optical type II sources). 

(5) Our limits on the Ha/H/3 ratio in 4 cases suggests that 
the optical emission in both narrow- and broad-line op- 
tically faint ULIRGs is heavily extincted (by at least 
2.4 magnitudes at the wavelength of Ha). Since the 
narrow-line region is also extincted, this suggests dust 
obscuration on large scales. 

(6) Given that these sources are bolometrically powerful, 
we would expect the AGN to be hosted in a large 
galaxy which is actively forming stars and optically 
bright. However, the UV continuum emission from 
the host galaxies is faint and in some cases may be 
extended. We suggest that these sources are examples 
of 'host obscured' AGN in which a large proportion of 
the obscuration is contributed on large (~ kpc) scales 
within the host galaxies. 



We intend to obtain similar observations of a larger sample 
of optically faint, luminous infrared galaxies to determine the 
fraction of type I to type II AGN more accurately and to get 
better constraints on dust extinction within these sources. 
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